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ABSTRACT 

An earlier study (Jones et al., 2000) compared NASA/NSO Spectromagneto- 
graph (SPM) data with spacecraft measurements of total solar irradiance (TSI) 
variations over a 1.5 year period in the declining phase of solar cycle 22. This 
paper extends the analysis to an eight-year period which also spans the rising 
and early maximum phases of cycle 23. The conclusions of the earlier work ap- 
pear to be robust: three factors (sunspots, strong unipolar regions, and strong 
mixed polarity regions) describe most of the variation in the SPM record, but 
only the first two are associated with TSI. Additionally, the residuals of a lin- 
ear multiple regression of TSI against SPM observations over the entire eight- 
year period show an unexplained, increasing, linear time variation with a rate of 
about 0.05 W yr~^. Separate regressions for the periods before and after 
1996 January 01 show no unexplained trends but differ substantially in regression 
parameters. This behavior may reflect a solar source of TSI variations beyond 
sunspots and faculae but more plausibly results from uncompensated non-solar 
effects in one or both of the TSI and SPM data sets. 


1. Introduction 

Accurate understanding of global solar variability is of clear astrophysical interest 
and is also vital to distinguish natural from anthropogenic causes of long-term changes in 
terrestrial climate. Daily averages of TSI observations from several spacecraft radiometers 
over the past two decades show clear rotational and solar-cycle variations (Pap and Frohlich, 
1999; Frohlich and Lean, 1998; Frohlich 2000). The modeling of this variability through 
comparison of the spacecraft measurements with spatially resolved solar observations from 
both ground- and space-based instruments is highly refined (Chapman et al., 1996; Foukal 
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and Lean, 1988; Fligge et al., 2000; Preminger et al., 2002). Two classes of solar features, 
dark sunspots and bright faculae, account for about 90% of the TSI variance. Although it 
is not yet clear whether the remaining discrepancies are observational or require additional 
sources of irradiance variability, the simple observation (de Toma et al., 1999) that TSI at 
the current solar maximum is very similar to the previous maximum while photospheric 
indicators of solar activity are lower in cycle 23 than in cycle 22 hints that more than solar 
activity is involved. This paper presents more detailed evidence to support this idea. 


2. Observations 

The results discussed here are based on two data sets. One of these is the set of daily 
full-disk magnetograms from the SPM which consist of strictly cospatial and cotemporal 
full-disk images of line-of-sight (LOS) magnetic flux, LOS velocity, continuum intensity, 
equivalent width, and central line depth derived from long-slit spectral polarimetry of the 
Fe I 868.8 nm line. The observations have 1.14" x 1.14" spatial pixels, approximately 42 
mA spectral pixels and a magnetic noise level of about 5 Gauss. Data for 2194 days from 
1992 November 21 - 2000 September 30 are reduced and analyzed here. A more complete 
description of the instrument and data is given by Jones et al. (2000), hereafter referenced 
as JBJW. 

The other data set is the composite measure of TSI compiled from various spacecraft 
from 1979 through the present by Frohlich and Lean (1998) and Frohlich (2000) who 
describe the individual observations and the techniques used to combine them into a 
single data stream. For this study we use version 23 of the composite TSI from the 
Physikalisch-Meterologisches Observatorium Davos World Radiation Center, Davos, 
Switzerland, which includes unpublished data from the Variability of Solar Irradiance 
and Gravity (VIRGO) experiment on the cooperative ESA/NASA Solar and Heliospheric 
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Observatory (SOHO). Data from both the composite TSI and the SPM were available for 
2126 days during the interval from 1992 November 21 to 2000 September 30. 


3. Analysis 

No adjustments or filtering were performed on the composite TSI data other than to 
select days for which both SPM and TSI observations were available. The reduction of the 
SPM data was carried out exactly as described in JBJW. Low-order, least-squares spatial 
fits of the daily data were applied to remove spurious instrumental and telluric effects 
as well as unwanted solar center-to-limb variations. Daily multi-dimensional histograms 
were computed from the five input images, and solar “features” were defined by restricted 
domains in the histogram variable space as described by JBJW. Area and contrast measures 
of each of the features were computed from the histograms as described below and in more 
detail in JBJW. Factor analysis using principal components extraction revealed the most 
important linear combinations of the original feature measures for explaining the variability 
of the SPM data, and a multiple regression of the TSI observations as a function of these 
factors was performed. This analysis procedure was meant to be exploratory, allowing 
efficient multiple passes through our data set to learn what important features in the SPM 
data are related to TSI without building in strong prejudices. We did indeed experiment 
with domain limits different from JBJW for the various features, but this did not affect the 
substantive results of our study. 

Specifically, we considered nine categories of solar features: weak-field; sunspot 
penumbra; sunspot umbra; bright (and dark) strong-field unipolar features in the central 
portion of the disk; bright (and dark) strong-field mixed polarity features in the central 
portion of the disk; bright (dark) strong-field features near the solar limb. Bright and 
dark features are distinguished by continuum (not bolometric) intensity contrast at 868.8 
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nm with reference to the fitted limb darkening curve. The defining parameter domains 
for these features both as used by JBJW and an alternative classification are given in 
Table 1. B is in Gauss, i and q are intensity and equivalent width contrasts with respect 
to local center-to-limb variation, and /x is the cosine of the heliocentric angle from disk 
center. “Unipolarity” , u, is defined here as the fraction of flux of one sign in a roughly 
supergranular area surrounding a given point on the disk and is divided into four bins 
in the multidimensional histograms. The revised feature classification was also used by 
Jones (1998) to search for possible thermodynamic differences between unipolar and mixed 
polarity features. The alternate and original classifications differ mostly in the division of 
unipolar and mixed-polarity features. Roughly speaking, strong unipolar regions correspond 
to active regions, including some of the enhanced network which trails from older active 
regions, while strong mixed polarity regions correspond to magnetic network outside of 
active regions. The revised feature definitions count more enhanced network as unipolar. 

As in JBJW, for each histogram feature we compute the fractional area (number 
of pixels relative to the total number of pixels) and differential contrasts in intensity, 
equivalent width, and central line depth relative to the weak-field pixels in the same domains 
of u and fx. As basis sets for subsequent factor analysis, we use two different subsets of 
these quantities. In both cases, we include fractional area and diflferential contrasts for 
each of the feature classes except weak fields where the differential contrasts are zero by 
definition. In the “original” basis set used by JBJW, we also include weak-field fractional 
areas for the whole disk, disk center unipolar and mixed polarity regions, and limb. In 
the “new” set, we include instead the weak-field fractional area and contrasts with respect 
to limb-darkening for intensity, equivalent width, and line depth for the whole disk. Our 
substantive conclusions do not depend on which feature definition or basis set is used. 

Altogether, 36 variables are computed for either basis set. Factor analysis, which 
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Table 1. Histogram Feature Domains 


Class 

1 B 1 

i 

Q 


u 

weak 

[0,8] 

all 

all 

[0,1] 

[0,1] 

penumbra 

(128,768] 

>-0.3 

>0.05 

[0,1] 

[0,1] 

umbra 

>128 

<-0.3 

>0.0 

[0,1] 

[0,1] 

disk center, original 

bright mixed polarity 

[8,256] 

>0.0 

all 

[V2/2 ,1.0] 

<0.75 

bright unipolar 

[8,256] 

>0.0 

all 

[y/2/2 ,1.0] 

>0.75 

dark mixed polarity 

[8,256] 

<0.0 

all 

[x/2/2 ,1.0] 

<0.75 

dark unipolar 

[8,256] 

<0.0 

all 

[x/2/2 ,1.0] 

>0.75 

disk center, revised 

bright mixed polarity 

[8,256] 

[0.0,0.1] 

[-0.1,0.05] 

[V2/2 ,1.0] 

<0.5 

bright unipolar 

[8,256] 

[0.0,0.1] 

[-0.1,0.05] 

[\/2/2 ,1.0] 

>0.5 

dark mixed polarity 

[8,256] 

[-0.1,0.0] 

[-0.1,0.05] 

[V2/2 ,1.0] 

<0.5 

dark unipolar 

[8,256] 

[-0.1, 0.0] 

[-0.1,0.05] 

[V2/2 ,1.0] 

>0.5 

bright limb 

[8,256] 

>0.0 

all 

[0,V2/2] 

[0,1] 

dark limb 

[8,256] 

<0.0 

all 

[0,V2/2] 

[0,1] 
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models the observed variables as a linear combination of unknown “common” factors plus 
unique variability (including noise), is used to reduce the dimensionality of the original 
basis by accounting for cross-correlations. Rather than iteratively solving the complete 
factor analysis model, we fix the unique variability at zero and identify the common factors 
as the principal components (eigenvectors) of the correlation matrix, select from inspection 
of the eigenvalue “scree” plot the factors with the six largest eigenvalues, and perform 
an orthogonal “varimax” rotation of these to maximize the column variance of the factor 
loading matrix while maintaining zero cross correlation between factors. The rotated 
factors, which are linear combinations of the original variable set, account for most of the 
variance in the SPM data and have relatively unambiguous interpretations in terms of solar 
features. A more detailed description of this process is given in JBJW. Finally, we perform 
a multiple regression with TSI as dependent and SPM factors as independent variables. 

Since the work of JBJW, the accuracy of the algorithm for computing the weak-field 
reference contrasts for intensity, equivalent width, and line-depth has been improved. In 
this paper, we report results from applying the improved algorithm to both the time 
period spanned by JBJW and that spanned by the complete data set. The factor loadings 
(correlations of the original variables with the factors) for the original and improved 
algorithms are highly correlated, and our conclusions would remain unchanged had we used 
the original algorithm. Although any combination of feature definitions and basis sets 
which we tried lead to the same conclusions, computations using the “new” basis variables 
extracted for the “revised” features show, by a small margin, the best combination of clear 
interpretability of the factor loadings and continuation of the factor structure over the 
entire time period and will be displayed in this paper. 
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4. Results 
4.1. Factor Analysis 

Figure 1 shows the fraction of explained variance for both the unrotated and rotated 
factors over the entire time period of the SPM observations. The six retained factors explain 
about 84% of the SPM variance, while the dominant three factors alone account for about 
69%. We assume that the remaining 30 components whose eigenvalues slowly decay to 
zero and which both individually and cumulatively account for a small fraction of the SPM 
variance are primarily indicative of noise in the data. Note that the orthogonal rotation of 
the components redistributes the explained variance mostly among the first three factors. 

Figure 2 shows the rotated factor loadings (i.e., the correlations of each factor with 
each of the original basis variables) in bar graph form for the entire observing period. 
Factors are numbered in decreasing order of explained variance. Over this period, the 
most important factor is highly correlated with fractional area and differential contrasts of 
intensity, equivalent width, and line depth for both bright and dark unipolar, disk-center 
features and bright limb features (with the exception of equivalent width for dark disk-center 
features) and fractional area of dark limb features and weak-field structures. This loading 
pattern represents spatial structures associated with active regions outside of sunspots 
together with enhanced unipolar network and at least roughly corresponds to faculae in the 
standard model of TSI variations. The second factor is associated with area and contrasts 
in sunspots, while the third factor depends on areas and contrasts for both bright and dark 
mixed-polarity, disk-center features (again with the exception of equivalent width for dark 
structures). One consequence of our revised feature definition is the cleaner separation of 
unipolar and mixed polarity regions in factors one and three than the original classification 
of JBJW while Factor 1, representing primarily faculae, shows stronger correlation with 
sunspots. As will be shown more quantitatively below, the dominant three factors closely 
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Fig. 1. — Proportion of explained variance as a function of factor number. Upper panel 
shows contributions for each factor individually, lower panel shows cumulative contribution. 
Open circles show the non-rotated factors (principal components) while plus signs denote 
the six rotated factors. 
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correspond to the factor structure described by JBJW for the limited time period. 

Although we do not yet understand how the loading patterns for the three minor 
factors relate to conventionally recognized solar features, they consistently appear with 
the improved analysis algorithm, new basis variable set, and revised feature definitions 
regardless of sampling period. Differential contrast of equivalent width for dark disk-center 
structures, regardless of unipolarity, is highly correlated with Factor 4; differential contrasts 
of intensity, equivalent width and line depth in dark limb features are related to Factor 
5; and Factor 6 is associated with contrasts of intensity (primarily) and equivalent width 
(secondarily) of weak magnetic elements. 

Figure 3 shows the rotated factor loadings for the original period treated by JBJW. 
Except for minor interchanges of order (Factors 1 and 3 as well as 5 and 6), Figures 2 and 3 
are visually very similar. The stability of the factor patterns is more quantitatively verified 
in Table 2 which shows the correlation matrix of the two factor patterns. The loading 
patterns are similar for all six of the retained factors and therefore are likely to refiect 
intrinsic properties rather than an accidental result of the sampling period. The different 
ordering over the two time periods refiects partly the varying importance of solar properties 
over different phases of the sunspot cycle and partly the near equality of explained variance 
among the relevant factors. 


4.2. Multiple Regression 

Table 3 shows the coefficients and their standard errors along with the fraction of 
explained variance (multiple R?) for three regressions of composite TSI as a function 
of the six rotated factors. The first covers the total time period 1992 November 21 - 
2000 September 30; the second and third regressions, following de Toma et al. (2000), 
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1992 Nov 21 - 2000 Sep 30 



Fig. 2. — Bar graph of the correlation coefficients between the rotated factors and original 
variables (factor loadings) for the entire SPM observation period. X-axis is labeled by original 
variable: i, q, d denote contrasts of intensity, equivalent width, and central line depth; n 
denotes fractional number of pixels (area); pen refers to sunspot penumbra, umb to umbra; 
unicb(d) refers to unipolar disk-center bright (dark) pixels; mixcb(d) refers to mixed polarity 
disk-center bright (dark) pixels; lb(d) refers to limb bright (dark) pixels; wk refers to weak- 
field pixels. 
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1992 Nov 21 - 1994 Mar 30 




Fig. 3. — Factor loadings for the initial SPM period treated by JBJW in the format of 2. 
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are independent fits of the periods before and after 1996 January 01, a division date 
corresponding roughly to the minimum separating cycles 22 and 23 as well as the beginning 
of observations from the VIRGO radiometers. The bottom row of Table 3 shows the slopes 
in W m~^ yr~^ of linear least-squares fits of the residuals (TSI - regression) for the three 
multiple regressions as a function of time together with the formal errors and fractions of 
the residual variances accounted for by the fits. The composite TSI (version 23) together 
with regression fits and residuals are plotted as functions of time in Figure 4. 

For the total period the six factors account for about 77% of the TSI variance, somewhat 
more than the similar analysis of JBJW but less than the best traditional analyses using 
PSI and PFI (e.g., Chapman et al., 1996). Moreover, as in JBJW, unipolar regions (faculae) 
and sunspots account for almost all the explained variance while strong-field mixed polarity 
(quiet network) features, although accounting for a significant fraction of the SPM variance, 
are not significantly correlated with TSI. For the declining phase of cycle 22, the second 
multiple regression accounts for noticeably less (71%) TSI variance while the third multiple 
regression accounts for nearly 89% of the TSI variance in the growth and maximum phase 
of cycle 23. For the latter two regressions, the major differences between TSI and SPM 
prediction occur on time scales of days to a rotational period, and the explained variance 
for cycle 23 is comparable to the best current fits obtained from other data. Although 
quiet network (Factor 3) has coefficients in the latter two regressions which are statistically 
different from zero at better than 3(7, the feature accounts for less than 1% of the TSI 
variance in either case. Although the present analysis is in qualitative agreement with 
JBJW, the use of improved analysis algorithms and particularly revised feature definitions 
(see Table 1) changes details; unipolar regions and sunspots account for comparable 
amounts of TSI variance in JBJW while in this analysis the regressions are dominated by 
unipolar features. 
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Table 2. Correlation of Factor Patterns 


Factors for 

1992 Nov 21 - 1994 Mar 30 

Factors for 1992 Nov 21 - 2000 Sep 30 
1 2 3 4 5 6 

1 

0.502 

0.253 

0.920 

-0.234 

-0.154 

0.173 

2 

0.494 

0.954 

0.135 

0.075 

-0.066 

0.038 

3 

0.914 

0.514 

0.024 

-0.115 

0.008 

0.059 

4 

-0.018 

0.107 

-0.521 

0.873 

-0.157 

-0.293 

5 

0.049 

-0.008 

0.059 

-0.308 

-0.388 

0.866 

6 

0.157 

0.044 

-0.021 

-0.146 

0.918 

-0.166 





Table 3. 

Multiple Regressions 




Factor 
(SPM B?) 

1992 Nov 21 - 
coeff error 

2000 Sep 30 
B? 

1992 Nov 21 - 
coeff error 

1995 Dec 31 
B? 

1996 Jan 01 - 
coeff error 

2000 Sep ; 
B? 

1 (0.30) 

0.346 

0.004 

0.685 

0.207 

0.007 

0.492 

0.376 

0.004 

0.768 

2 (0.20) 

-0.112 

0.004 

0.073 

-0.188 

0.008 

0.204 

-0.136 

0.004 

0.101 

3 (0.18) 

-0.010 

0.004 

0.001 

-0.022 

0.006 

0.006 

-0.037 

0.004 

0.005 

4 (0.06) 

0.035 

0.004 

0.007 

0.012 

0.006 

0.002 

0.047 

0.004 

0.012 

5 (0.05) 

0.020 

0.004 

0.002 

0.003 

0.005 

0.000 

0.008 

0.005 

0.000 

6 (0.04) 

0.009 

0.004 

0.000 

0.003 

0.004 

0.000 

-0.035 

0.007 

0.002 

Total (0.83) 



0.768 



0.705 



0.889 

Trend 

0.049 

0.002 

0.292 

-0.005 

0.006 

0.001 

0.006 

0.003 

0.003 
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Fig. 4. — TSI (top), multiple regression fits (middle), and residuals (TSI -fits; bottom) as 
a function of time in years from 1992 November 21 - 2000 September 30. All abscissae 
are in W m“^. Heavy line in upper curve of lower panel shows linear least-squares fit to 
the residuals for the entire period. Lower curves of middle and bottom panels show fits 
and residuals for regressions before (after) 1996 January 01, displaced by 1 W m“^, as solid 
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A new feature appears in the current regression analysis for the total time period 
which can be seen both in the bottom row of Table 3 and in the bottom panel of Figure 
4. In addition to noise on short time scales (days- weeks) , the plot shows systematic linear 
variation with a slope of about 0.05 W m“^ yr“^ The fit accounts for 29% of the variance 
in the residual curve, i.e., about 7 % of the TSI variance, suggesting that if this trend were 
removed from the TSI data, multiple would increase to about 0.84. Over an 8-year 
period this trend is a significant fraction of the TSI variation from solar minimum to 
maximum. One can also see from Table 3 and Figure 3 that there are no significant residual 
trends if the fit is segmented at 1996 January 01. However, the regression coefficients of the 
dominant factors for the two time periods differ from each other well beyond their formal 
errors in the sense that the composite TSI after 01 January 96 depends more strongly on 
Factor 1 (faculae) and less strongly on Factor 2 (sunspots). As discussed below and in 
agreement with de Toma et al. (2000), these results imply either that there are systematic 
observational errors in one or both the TSI or SPM data sets or that an additional solar 
source of irradiance variation exists which is not detected in the SPM data. 


5. Discussion 

The two major conclusions of JBJW are shown to extend to a decadal time scale. 
First, unipolar magnetic areas associated with active-region and active-network faculae 
and sunspots dominate the correlation of SPM observations with TSI. Second, strong-field, 
mixed polarity regions (quiet network), although contributing substantially to the total 
variance of the SPM record, are effectively uncorrelated with TSI. The duration of our 
observations is not extensive enough to exclude quiet network as an important long-term 
source of solar irradiance variations. However, our data does span most phases of the 
sunspot cycle, and the correlation of TSI with mixed-polarity network is low enough to 
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imply that any such contribution is likely to occur only on temporal scales considerably 
exceeding a solar cycle. Our analysis does not provide further insight into why TSI is poorly 
correlated with mixed polarity network. We note that Jones (1998) was unable to find 
any indications of thermodynamic differences between unipolar and mixed polarity regions 
with comparable LOS fluxes and speculate, as did JBJW, that quiet network is associated 
with “magnetic carpet” (Title and Schrijver, 1998) which disappears and renews over times 
much shorter than a solar rotation period. 

The linear temporal trend in the residuals of the multiple regression over the entire 
interval spanned by our data is consistent with a non-magnetic solar source of TSI variation 
or with long-term, systematic, instrumental effects. We are unaware of any likely source of 
long-term instrumental trends in the SPM measurements of LOS field but plan to check 
on this possibility by comparing SPM observations with other measurements. Even if such 
a trend were found, however, we note that magnetic fiux enters into our analysis only 
through the division of histograms into features and is not otherwise used in determining 
the measures which comprise the basis set for the factor analysis and multiple regression. 
Thus, a change in SPM sensitivity over time would modify the temporal variation of each 
factor but, to the extent that our representation is complete, would not affect the overall 
regression. There may, of course, be irradiance changes associated with magnetic features 
below the sensitivity or resolution limits of the SPM which would be undetected in our 
measurements. 

In any case, a trend produced by a continuous phenomenon, such as might be expected 
from the changing sensitivity of a single instrument or a long-term solar variation, should 
be observed in regressions of subsamples of the data. However, evidence for long-term 
systematic variation of the residuals in the segmented regressions is absent which suggests 
that the trend for the entire time period may be an artifact of forcing a continuous fit 
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to disjoint intervals. Although the fit for either interval would be consistent with the 
two- component activity model of TSI variation, the dominant regression parameters for 
cycles 22 and 23 are distinctly different. It is possible to construct scenarios for which the 
response of actual TSI to solar activity, as measured by SPM, is different in cycle 23 than it 
was in cycle 22. For example, Livingston (2002, private communication) finds that sunspot 
umbrae in the 1.6 micron lines near the maximum of cycle 23 show weaker magnetic 
fields and intensity contrasts than do those drawn from a comparable sample during the 
maximum of cycle 22. These differences might be undetected in SPM measurements due 
to substantial photospheric stray light and would appear as a weaker sensitivity of TSI to 
observed sunspot properties in cycle 23. Unfortunately, it is not possible to test this idea 
directly since the SPM was not operational during the maximum of cycle 22. Similarly, if, 
in cycle 23, unipolar network and faculae were brighter in the UV for given photospheric 
properties, then an enhanced apparent response of TSI would be observed. Such a scenario 
would be consistent with the work of Unruh et al. (1999) and Preminger et al. (2002) 
who suggest that lines are responsible for the bulk of the cyclic variability of TSI as well 
as with observations of the Mg II h and k resonance lines which tend to show comparable 
levels of emission in cycles 22 and 23. However, a source of chromospheric heating which is 
independent of photospheric magnetic properties detectable by the SPM would be required. 

A simpler and perhaps more likely explanation involves uncompensated systematic 
instrumental effects. While there are no documented changes in SPM instrumentation or 
sensitivity between cycles 22 and 23, the observed composite TSI is sensitive to changes in 
the ensemble of operational radiometers, and one such change occurred on 1996 January 18 
with the onset of VIRGO observations. Thus, the composite TSI, even though it has been 
carefully prepared to compensate for individual instrument sensitivities, plausibly records 
TSI with a different “gain” in cycle 23 than in cycle 22. Note that the SPM data fit the 
cycle 23 (predominantly VIRGO) observations best. 
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6. Future Work 

Our histogram-based technique tends to account for less TSI variance than the best 
representations based on other ground-based data, particularly in cycle 22. To understand 
better what aspects of our analysis lead to this result, we plan to compare the spatial 
structures which represent the factor patterns outlined in section 4.1 with other methods 
of feature classification, particularly those developed by Harvey and White (1999), Turmon 
et al. (2002), and Preminger et al. (2000). The work of Turmon et al. (2002) in particular 
may suggest better posed multidimensional alternatives to our histogram classification. 
Functional relationships between all the SPM observables derived from histogram analysis 
can be cross-checked with predictions from extant models such as those constructed by 
Fontenla et al., (1999), Solanki and Unruh (1998), Fligge et al. (1998), and Unruh et al., 
(1999) to explore physical differences between features. 

We have not attempted to separate quantitatively rotational, decadal, and other 
possible temporal scales in the representation of either TSI or SPM data. However, we plan 
to apply singular spectrum analysis (Pap and Varadi, 1996) and possibly other techniques 
such as wavelet analysis in the near future. This may help to understand what parts of the 
unexplained variance are attributable to different scales of variation and to clarify the lack 
of correlation between mixed polarity features and TSI. 

Finally, longevity and continuity are essential ingredients for any observational study 
of solar irradiance. Within calendar year 2003, the SPM will be retired from service 
to be replaced by the Vector Spectromagnetograph (VSM), a part of NSO’s Synoptic 
Optical Long-term Investigations of the SUN (SOLIS) instrumental package. The VSM 
will continue the observational repertoire of the SPM, albeit using a different spectral line, 
and will provide hitherto unavailable information regarding the vector magnetic field over 
the entire solar disk. Thus, important parts of our future work will be to complete analysis 
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for the SPM data set, to ensure that the transition between the two instruments produces 
minimal discontinuity in the synoptic magnetogram record, and to develop methods for 
including the more complete physical description of the solar atmosphere provided by the 
VSM into our irradiance analysis. 

The authors acknowledge useful discussions with G. Chapman, K. Harvey, W. 
Livingston, J. Pap, M. Turmon, and S. Walton as well as use of computer facilities at the 
University of Arizona for carrying out the factor and multiple regression analyses. This 
research was partially supported by NASA Supporting Research and Technology tasks 
344-12-52-14 and 344-12-52-19. NSO/KPVT data used here were produced cooperatively 
by AURA/NSO, NASA/GSFC, and NOAA/SEC. 
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